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A special model of a massless vector-field is presented, which has an extra modified- 
gravity-type interaction term in the action. The cosmology of the model is studied with 
standard noninteracting relativistic matter added. It is found that this cosmology can 
have an early phase where the vector-field starts to compensate a (Planck-scale) cos- 
mological constant and a late Friedmann-Robertson— Walker (FRW) phase where the 
relativistic-matter energy density dominates the dynamic vacuum energy density. 
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1. Introduction 

The cosmological constant problenP^l ig perhaps the most important outstand- 
ing question of modern physics. In a nutshell, the problem is to explain why 
the energy scale corresponding to the measured value of the cosmological con- 
stant A is negligible compared to the energy scales of elementary particle physics, 

|A|l/4 « £;qcd « -Bcw « -Bpianck ^ lO^^ GcV. 

The cosmological constant A can be canceled dynamically and without fine- 
tuning in a particular massless- vector-field modelP^l g^t this model runs into two 
obstacles. First, the local Newtonian dynamics is ruined.^ Second, the final phase 
of the model universe expands so rapidly (scale factor a on t) that any standard 
matter contribution {pM oc 1 /a^ with p — 4 for relativistic matter or p — 5 for non- 
relativistic matter) becomes irrelevant compared to the decreasing remnant vacuum 
energy density {pv oc 

Following previous work on the g-theory 

approaclJHmHMinillini 

to the main 

cosmological constant problem, an extended model with two massless vector-fields 
has been found ,^1^^^ which evades the first obstacle regarding the local Newtonian 
dynamics. Here, we present another model with a single vector-field and a modified- 
gravity-type interaction term, which circumvents the second obstacle by having a 
final phase with slower expansion (a oc t^^^). 



2. Model and Ansatz 



Consider the following model of a massless vector- field Aa (x) with an effective action 

R 



So 



-9 



+ e(Q3) + A + /(i?) A^A'^ + Cm 



(la) 



16nG 

where R is the Ricci scalar, A the effective cosmological constant, and f{R) a 
function to be specified later. For the moment, the contribution of other matter 
fields is set to zero. Cm — 0. The g-theory-type scalaj™ entering the above 
action is defined as follows: 



Ux)^g''^'{x)A^,p{x), 



(lb) 



with the spacetime coordinate (x) — {x^, x^, x^, i) as argument and the semicolon 
standing for covariant derivation. The particular contraction (jlbp was mentioned in 
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Ref. m but not studied in detail. 

In the following, we consider a spatially-flat Friedmann-Robertson-Walker 
(FRW) universe with scale factor a{t) for cosmic time t and take the isotropic 
AnsatM^ for the vector-field: 

(g^pit)) = (diag[l, ~a^{t), ~a\t), -a\t)]'^ , (2a) 

A^{t)^Ao{t)6l. (2b) 
The Ansdtze ^ give 

R^~6{H + 2H^), (3a) 

Q3= Ao + 3HAo, (3b) 

where the overdot stands for differentiation with respect to the cosmic time t and 
the Hubble parameter H{t) is defined by H{t) = a{t)/a{t). 

Note that the variable (Qs)^ = {Aq + 3H Aq)^ differs from the variable (Qi)^ = 
Aa-p A"-^ = Al + 3H'^ Al discussed in Ref s. 151 and [HI For the statics of g-theory, 
the precise realization of the g-type vacuum variable is irrelevant [here, either 
or Qi]. But for the dynamics, the realization matters [here, giving a FRW universe 
with tH{t) approaching either the value 1/2 or the value 1]. See App. A of Ref. [T^ 
for a brief introduction to the basic ideas of g-theory. This introduction to g-theory 
may give the reader some background information but is not needed to understand 
the present article, which is entirely self-contained. 

3. Reduced Field Equations 
3.1. Vector-field equation 

The variational principle applied to the vector-field of action ([Ta|) gives the following 



-2/(i?)A„. (4) 



equation: 

The Ansdtze ([2]) reduce this field equation to a single ordinary differential equation 
(ODE), 

CAo + {3HC + OAo+[3HC + 3HC- .f{R)]Ao = 0, (5) 
with C= {2 Qsr^e/dQs. 
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3.2. Generalized FRW equations 

The energy-momentum tensor of the vector-field Aa follows from the variation of 
the action with respect to the metric field ga[j{x). The result is 



T. 



de 
dQ~3 



\dQs 



gap + 2 V^) 



de 
dQ3 



-2 (^A^ f R^p -^A^g^p + f A^Ap - V„ (/' A') + g^p (/' A^)^ , (6) 

where the round brackets around spacetime indices denote symmetrization and 
the prime on / stands for differentiation with respect to R. With the vector-field 
equation ([5]), the energy-momentum tensor ^ becomes 



al3 



e-Q3 



de 
dQa 



-2 (^A^ f R^p + lA^fg^p-f A^Ap - V, (/' A^) + g^p (/' A^)^ . (7) 

From the gravitational field equations, the Ansdtze ([2]), and the energy- 
momentum tensor ([7]), the generalized FRW equations are 



2H + = (Stt G) [a - P{A) - wm Pm] , 



(8a) 
(8b) 



having added a further matter contribution from a nonzero term Cm in the ac- 
tion Pa|) . with W]^J = Pm/pm the constant equation-of-state (EOS) parameter of 
the corresponding homogeneous perfect fluid. The vector-field contributions on the 
right-hand sides of ([5]) are given by 



piA) 
PiA) 



de 
do's 



3iH + H') f Al + lf AI-3H ^^[r Al 



, (9a) 



dQs 



-2 



{H + 3H^)fAl--fAl 



(9b) 
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4. Special Vector-Field Model 

4.1. Linear f{R) and quadratic €{Qs) 

The special model considered in this article has a linear function f{R) in the inter- 



action term of the action ([Ta|) . 

J{R) = kR, (10a) 

and a quadratic vector-field function (.(Qs), 

e(Q3) = CoQl, (10b) 

for nonzero dimensionless coefficients n and ^o- In order to be able to cancel a 
nonzero cosmological constant, these two coefficients must obey the condition 

sgn(6 K + 5Co) = sgn(A) , (10c) 

in terms of the sign function sgn(a;) = x/\x\ for x and sgn(O) = 0. The particular 
form of condition (|10c|) will be derived at the end of this subsectionQ 

The Ansatz ()10a|) reduces the energy density and the isotropic pressure from ^ 
to the following expressions: 

p{A) = +e-Q^^^-Qn{H^Al + Hj^ Al^ , (11a) 



Similarly, the Ansatz (|10b[) reduces the ODE (O to the following equation: 



3ij-l/(i?) 
so 



Ao = 0. (12) 



Inserting the function (|10ap in the last ODE and using ((3a|) gives 

io + 3iJio + 3ijylo + 6-^ (i/ + 2i72) = 0. (13) 

Co 

Equation (fT5|) is the core result of this article, as will become clear in Sec. 14.21 
For t — > cx), the asymptotic solution of ()13p is given by 

Aasymp(t) = , ifasymp(i) = (1/2) , (14) 



'^Condition lllOcl l is purely technical and appears because of the simple e-function used in l llObl l. 
For an appropriate, more complicated, e-function (cf. Refs. [6l rfl l8l [T3t . the cancelation of the 
cosmological constant A can be expected to hold for any sign of A. 
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with a constant to be determined by the solution of the combined field equations 
(see Sec. 14. 3p . Substituting this asymptotic solution into the energy density p{A) 
and the pressure P{A) from ([TT|) . we obtain 



PA(^as.vmp) = +A + 



-PA(^asymp) — ^ ' 



de 
de 



Q3=5 0/2 2 



15 —2 

K0 , (15a) 



Q3=5 0/2 2 



^^/^6l^ (15b) 



having added the contributions from the cosmological constant A in the action (|la[l . 
Observe that pA(^asymp)+^A(^asymp) — 0. In fact, (fTS)) has the particular structure 
of g-theory.fSEIHl 

It is even possible to define a new effective e-function by making 

2 

a Qa-independent shift, IcsiQa) — f^iQa) ~ (15/2) k 9 . 

Equation psp also explains condition (|10c|) : for the simple quadratic function 
(jlObp , it is possible to cancel A with an appropriate value of the linear 6'-coefficient 
from p4)) only if condition (|10cl) is satisfied. As mentioned before, the actual value 
6 follows from the solution of the combined field equations. Remark, finally, that 
the expression on the right-hand side of (|15ap may be called the 'dressed' cosmo- 
logical constant if A from the action ([Ta|) is considered to be the 'bare' cosmological 
constant. 



4.2. Heuristic argument 

Returning to the vector-field ODE p^ . we can give the following heuristic argument 
for the appearance of the FRW-like asymptotic solution , even with a nonzero 
cosmological constant A present in the action (|la|) . The physically relevant case has 
both coefficients k and (o nonvanishing, but, for completeness, also the other cases 
will be briefly mentioned. 

For K 7^ 0, Co 7^ 0, and a linear time-dependence of the vector component Aolt) cx 
t, the H^Aq term in (IT5|) excludes having H{t) = const. Purely for dimensional 
reasons, assume H{t) = j t^^ with a numerical constant 7 > [excluding the 
case of a static or contracting Universe with 7 < ] . The first three terms of the 
left-hand side of ([T3| then cancel by themselves. This leaves the k term in (|13p . 
which vanishes for 7 = 1/2 corresponding to a radiation-dominated FRW universe. 
Indeed, the functions pi)) allow for a compensation of the cosmological constant A, 
as shown by ([T5|) and the final ODEs to be presented in Sec. 14.31 
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For K = and Co 7^ 0, the ODE (IT51) also has the de-Sitter solution with constant 
Aq and H. As said before, precisely this solution is forbidden by having the k H^Aq 
term in ([13]) if k 7^ 0. 

For K ^ and Co = 0, implies Aq = (i.e., 6* = 0) for generic H{t) and the 
cosmological constant can no longer be canceled, as exemplified by (I15|) . 

For K = Co = in the Ansdtze ([TO)) . the vector-fields have disappeared from the 
action (jlap altogether and the cosmological constant problem remains unsolved. 

Returning to the case of k 7^ and Co 7^ 0, the above heuristic argument for 
obtaining a radiation-dominated FRW universe, thus, relies on having an action 
employing the Qs field [which gives the 3i7 Ao + 3ijv4o combination in ([T^] and 
the Ricci scalar (which vanishes precisely for H = ^t^^). 

An entirely different question is whether or not solution (jl4p appears dynami- 
cally as an attractor. Here, this question will be addressed numerically, leaving a 
proper mathematical treatment for the future (see Note Added). 

4.3. Additional relativistic matter and dimensionless ODEs 

From now on, we also consider an additional standard-matter component with en- 
ergy density p]\j and pressure Pa/. These and other dimensional variables can be 
replaced by the following dimensionless variables: 



if appropriate powers of the (Planck-type) energy scale (Stt G)"^/^ are used, without 
further numerical factors. Henceforth, an overdot stands for differentiation with 
respect to r, for example, /i(t) = a{T)/a{T). 

The previous results (jSbp and (|13p . together with the extra standard-matter 
component, give the following dimensionless ODEs: 



{A, e, pm, Pm] -> {A, e, tm, Pm] , 



(16a) 



{t, H, Qs, Aq} {r, h, qs, v} , 



(16b) 




(17a) 



v + 3hv+ [ih + 6{K/Ca){h + 2h^)] u = 0, 



(17b) 



fj\/ -I- 3 (1 4- wm) hrM = 0, 



(17c) 
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where the last equation describes the evolution of the additional homogeneous 
perfect fluid of the standard-matter component, with constant EOS parameter 
wm = PmI^m ~ 1/3 for ultrarelativistic particles. The corresponding generalized 
Friedmann equation from (jSap is given by 



Observe that the right-hand side of , evaluated for the asymptotic solution (|14p , 
has canceling 0(t°) terms for one particular value of , provided condition (jlOc|) 
holds. Specifically, the two possible 6 values are given by 



in terms of the dimensionless cosmological constant A = (Stt G)^ A. The actual sign 
of 9 will be determined by the initial boundary conditions. 

In fact, the ODEs ((T7)) are to be solved with boundary conditions on v{t), v{t), 
h{T), and ^^/(t) at t = Tstart, where these particular function values must satisfy 
the constraint equation (jl8p . The corresponding physical quantities must be small 
enough for classical gravity to be relevant, e.g., rA/(Tstart) ^ 1- 

4.4. Numerical results 

Numerical calculations have been performed for the model defined by ([T]) and (jlOp , 
with the dimensionless cosmological constant A = 0.02 and model parameters k = 
—Co/2 — —1/2. The results are presented in two figures, the first without and the 
second with dynamical effects from the vector-field. 

With special boundary conditions to ensure the exact vanishing of the vector- 
field. Fig. [1] shows that the cosmological constant A rapidly dominates the matter 
component and that exponential expansion sets in. Asymptotically, this model uni- 
verse approaches de Sitter space having h — tm = 0. The behavior of ft.(r) and 
^A/(''") shown in Fig. [1] is completely standard and can also be obtained analytically. 

With boundary conditions of the vector-field in an appropriate domain. Fig. [5] 
shows that the vector-field is able to compensate a Planck-scale cosmological con- 
stant A = 0(1) and that Minkowski spacetime is approached asymptotically, 
/i(t) — > for r — )• oo. The numerical solution v{t) asymptotically takes the form 
(fT4| , with the precise linear coefficient needed for the complete compensation of the 



3h^ = X-Ca{v + 3hvY + tm - 6 k {h^ + 2 hv v) . 



(18) 




(19) 
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Fig. 1. Numerical solution of ODEs 1)17^ for parameters A = 0.02, = Ij ^ = —1/2, and 
= 1/3. The boundary conditions are i>(10) = i)(10) = 0, tmCIO) = 0.02, and h{10) = 0.115470. 
With these boundary conditions, the ODEs have the exact solution v{t) = for r > 10. 




Fig. 2. Same as Fig. [T] but now with nonzero starting values for the vector-field, v{lO) = 1 it 
0.1 and i)(10) = 0, where the dashed curves correspond to the starting value d(10) = 0.9. The 
corresponding starting values for h follow from ljl8|l and are, respectively, h{10) = 0.0624391 and 
h(10) = 0.0713376. The scaling in the ti-panel uses the constant from l|19| l. which, for A = 0.02, 
takes the value 9 = \/2 A/5 = 1/(5 \/5)- The numerical solutions of v{t) shown in the left panel 
approach asymptotically the function Bnum t, with Snum = 9 within the numerical accuracy of the 
calculation (the rescaled functions plotted in the left panel run towards the value 1 as r — > oo). 
The corresponding numerical solutions of /i(t) shown in the middle panel approach asymptotically 
the function 7numT^^, with 7num = 1/2 within the numerical accuracy of the calculation (the 
rescaled functions plotted in the middle panel run towards the value 1 as r — > oo). 

cosmological constant A. It needs to be emphasized that this hnear coefhcient of 
v{t) is not put in by hand but arises dynamicahy; cf. Ref. |8] See also the remark 
below (fT8|) and the details given in the caption of Fig. [2] 

The final phase shown in Fig. [2] corresponds to an FRW universe {h ^ ^ ■'"~^) 
dominated by relativistic matter. The solid curve in the right panel of Fig. [2] cor- 
responds to an asymptotic universe with vm ^ 1.09 and ry = (3/i^ — r^f) ''^ 
— 0.34r~^ (hence, ratio rjvj/l^'yl ^ 3) and the dashed curve to an asymptotic uni- 
verse with tm ~ 0.687 and ry ^ 0.063 (hence, ratio VM/ry ^ H)- For 
these numerical solutions, the action density term RAq cx —6 {h + 2 h"^) has also 
been found to drop to zero faster than t~^. 
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5. Discussion 

The vector-field model defined by ([1]) and (|T0)) provides an attractor solution which 
compensates an arbitrary positive cosmological constant A and gives a final uni- 
verse with Hubble parameter Hit) — \t~^ . This final state resembles a standard 
radiation-dominated FRW universe^ 

If the model is extended by the addition of a standard noninteracting-relativistic- 
matter component, the final state can be a genuine radiation-dominated FRW uni- 
verse with a subleading (time-dependent) vacuum-energy-density component. Fig- 
ure [2] shows two possible model universes, which start out with equal matter and 
vacuum energy densities, pM = A ~ (-Epianck)^, but end up with the matter com- 
ponent dominating, pu > \pv\, both components pM and \pv\ decreasing as 
(specific numbers are given in the last paragraph of Sec. 14. 4p . 

In this way, there is a more or less realistic physical description of the earliest 
cosmological phase as a radiation-dominated FRW universe with a dynamically 
canceled cosmological constant. Possibly, this description needs to be augmented 
with the effects from inflation.^l^l 

As discussed in Refs.l^ lTUlfTTl quantum-dissipative processes can be expected to 
lead to a further (exponential) reduction of 1/31/(^)1 in the very early universe (kT 3> 
Ecw ^ TeV) and related processes at the electroweak scale can perhaps generate 
a finite remnant value of order pv{oo) ^ (i?f^/i?pianck)^ ^ (meV)^. Alternative 
explanations of the remnant vacuum energy density in the g-theory framework 
have been reviewed in Ref. [121 

To conclude, both obstacles mentioned in the Introduction have been dealt with 
separately, the first in our previous article of Ref. [13] and the second in the present 
article. It remains to find a joint solution, provided such a solution exists. 



''Most likely, the model can be extended to allow for a similar compensation of a cosmological 
constant A of arbitrary signl«|7|«|13|For the case of positive A, also a final cosmological phase with 
Hit) = (2/3) can be obtained by replacing eiQs) in ^ with e(Qi, Q2) = Ci (Qi)^ + C2 iQl f, 
for the contractions (Qi)^ = A^.^ A"''^ and (Q2)^ = ^a;/3^'^'"i E^nd by taking the coefficient 
K = — (Ci + C2)/2 in {TOaJ- This type of model may be relevant to inflationary effects in vector- 
field theories with a dynamically canceled cosmological const ant At this point, it should be 
mentioned that certain scalar-tensor theories (with a fundamental scalar field 0) have also been 
argued to give both a compensation of the cosmological constant and a final FRW-like phase.fIS] 
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Note Added 

For the case of A > and k = — Co/2 = —1/2, it is possible to prove the existence 
of an asymptotically stable (attractor) solution. The proof relies on an appropriate 
change of variables, knowledge of the solution rM{a), and the Poincare-Lyapunov 
theorem [Theorem 7.1 in Ref. [TB]. Most likely, the attractor solution can be estab- 
lished rigorously also for general k ^ and Co 7^ with 6 k + 5 Co > 0. 

The present article is the second of a trilogy of articles, the first one being 
Ref. [131 The third article of the trilogy (Ref. 117]) finds the joint solution mentioned 
in the last sentence of Sec. [5j This third article also gives a detailed mathematical 
discussion of the attractor behavior in the type of vector-field models considered. 
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